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Abstract 

A  new  fluid  model  including  electron  pressure  is 
described  for  the  electron  flow  layer  in  a  magnetically 
insulated  transmission  line.  Non-zero  values  of  the 
electric  field  at  the  cathode  are  allowed  so  that  the  model 
can  treat  both  emission  and  re-trapping  of  flow  electrons. 
The  model  recovers  previous  results  when  the  electric 
field  at  the  cathode  is  taken  to  zero.  The  resulting 
generalized  MITL  flow  equations  are  discussed  and  initial 
solutions  are  presented. 


I.  INTRODUCTION 

Many  modern  pulsed  power  generators  use 
magnetically  insulated  transmission  lines  (MITL)  to 
couple  power  between  the  driver  and  the  load  [1],  In  an 
MITL  the  electric  field  stress  on  the  cathode  exceeds  the 
vacuum  explosive-emission  threshold  and  electron 
emission  occurs.  For  sufficiently  high  current,  emitted 
electrons  are  magnetically  insulated  from  crossing  the 
anode-cathode  gap  and  flow  axially  downstream  in  the 
direction  of  power  flow.  The  return  current  from  the  total 
anode  current  Ia  is  divided  between  current  Ic  flowing  in 
the  cathode  and  current  If  flowing  in  the  vacuum  electron 
layer,  i.e.,  If  =  Ia  -  Ic.  As  a  result  of  the  electron  flow  in 
vacuum  between  the  electrodes,  the  impedance  of  the 
MITL  is  altered  and,  thus,  the  power  flow  coupling 
between  the  machine  and  the  load  changes.  The  effective 
impedance  is  best  described  by  the  flow  impedance  Zf 
[2,3],  In  a  dynamic  system  where  the  voltage  and 
currents  are  changing  in  time,  Zf  also  varies  [4]. 

The  objective  of  the  work  is  to  efficiently  and 
accurately  simulate  power  flow  in  systems  with  a  MITL 
using  a  simple  TLC  rather  than  a  more  computer  intensive 
particle-in-cell  code.  A  generalized  model  for  flow 
impedance  is  developed  for  incorporation  into  a 
transmission  line  code  (TLC),  such  as  Bertha  [5].  The 
model  describes  both  self-limited  flow  as  the  pulse 
initially  propagates  down  the  MITL  toward  the  load,  as 
well  as,  the  subsequent  electron  power  flow  along  the 
MITL  after  the  pulse  encounters  the  load.  Additionally, 
for  low  impedance  loads  (i.e.,  with  impedance  less  than 
the  self-limited  impedance  of  the  MITL),  the  flow  is 
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Figure  1.  Schematic  of  MITL  flow  in  planar  geometry. 


modified  by  the  wave  reflection  off  the  load  and  the 
percentage  of  the  return  current  in  vacuum  electron  flow 
decreases.  This  phenomena  is  known  as  electron  re¬ 
trapping  [6]  and  must  be  included  in  a  complete 
description  of  MITL  flow. 

The  approach  is  to  develop  a  new  fluid  model  for  the 
electron  layer  which  includes  electron  pressure.  This 
modification  extends  MITL  flow  theory,  allowing 
solutions  with  non-zero  electric  field  at  the  cathode  (i.e., 
Ec  ^  0)  so  that  electron  emission  and  re-trapping  can  be 
treated.  A  well-defined  flow  impedance  is  also  derived. 
In  addition,  the  pressure  term  is  chosen  appropriately  to 
recover  previous  Ec  =  0  solutions  so  that  the  results  of  the 
new  rescaled  MITL  theory  [3]  can  be  applied.  This  newly 
extended  MITL  flow  model  may  then  be  combined  with 
time-dependent  field  equations  [4]  to  build  a  dynamic 
model  for  MITL  flow  in  a  TLC.  However,  to  complete 
the  TLC  model,  a  robust  numerical  technique  must  be 
constructed  for  solving  these  new  MITL  flow  equations 
and  techniques  must  be  developed  to  treat  special 
situations  such  as  impedance  transitions,  adders,  non¬ 
emitting  regions  (knobs),  and  load  coupling,. 


II.  GENERALIZED  MITL  FLOW  MODEL 

In  Fig.  1,  negative-polarity  MITL  flow  is  illustrated 
where  ra  and  rc  would  represent  the  anode  and  cathode 
radii  in  cylindrical  geometry.  For  simplicity,  the  model  is 
developed  in  planar  geometry,  considering  a  transmission 
line  element  of  length  £,  width  w,  and  anode-cathode  gap 

d.  The  three  basic  equations  for  the  fluid  model  are  the 
radial  momentum  transport  equation,  Poisson's  equation, 
and  Ampere's  law,  which  are  given  by 
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Here,  T  is  the  xx  component  of  the  electron  pressure 
tensor,  p  is  the  electron  charge  density,  vz  is  the  axial 
electron  fluid  velocity,  E  and  B  are  the  electric  and 
magnetic  fields,  and  e0  and  |U0  are  the  permittivity  and 
permeability  of  free  space,  respectively.  As  expressed  in 
Eqs.  (1)  -  (3),  p,  vz,  T,  and  B  are  all  positive-definite 
quantities,  while  E  is  positive  at  the  anode  but  can  be 
positive  or  negative  at  the  cathode.  Note  that,  because  the 
time  derivatives  have  been  assumed  negligible,  this  is  a 
quasi-equilibrium  model.  The  subscripts  c,  5,  and  a  will 
denote  quantities  evaluated  on  the  cathode  at  x  =  0, 
electron  layer  edge  at  x  =  ds,  and  anode  at  x  =  d, 
respectively.  The  vacuum  impedance  of  the  line  is  given 
by  Z0  =  c|i0d/w  where  c  =  l/(e0|U0)lc  is  the  speed  of  light. 
Also,  anode  and  cathode  fields  are  given  by  E  c  =  Z0cQa./d 
and  Bac  =  Z0Iac/ed.  where  Qac  and  Iac  are  the  charge  per 
unit  length  and  current  (on  the  anode  and  cathode), 
respectively. 

The  voltage  in  the  gap  is  given  by 


V(x)  =  jjE(x')dx'  ,  (4) 

where  Vc  =  V(0)  =  0  and  the  full  voltage  across  the  gap  is 
V  =  V(d)  (with  the  subscript  a  suppressed  for  simplicity). 
Similarly,  the  magnetic  vector  potential  in  the  gap  is 
given  by 


A(x)  =  j*B(x')dx'  ,  (5) 

where  Ac  =  A(0)  =  0  and  the  vector  potential  at  the  anode 
is  A  =  A(d)  (again  with  the  subscript  a  suppressed).  Both 
V  and  A  will  be  needed  for  a  complete  description  of  the 
problem. 

Because  electron  space-charge  and  current  are 
distributed  in  the  line,  there  is  not,  in  general,  a  well- 
defined  impedance.  The  electric  flow  impedance  Zf  and 
the  magnetic  flow  impedance  Zm  are  defined  by  the 
centroid  of  the  space-charge  and  current,  respectively  [2]. 
By  assuming  here  that  both  the  charge  density  and  axial 
fluid  velocity  are  uniform  in  the  flow  layer,  a  well-defined 
flow  impedance  is  derived  with  Zf  =  Zm.  This  assumption 
also  makes  integration  of  the  above  equations 
straightforward.  However,  without  the  3T/3x  pressure 
term,  Eq.  (1)  cannot  be  satisfied  under  this  assumption. 
Unlike  the  parapotential  flow  model  [7]  in  which  all 
electrons  move  axially  in  straight-line  orbits,  this  model 
permits  various  orbit  types  [8],  The  electron  pressure  is  a 
consequence  of  the  motion  perpendicular  to  the  axial  flow 
which  results  from  the  distribution  of  electrons  with  these 
various  orbits. 

The  MITL  equations  derived  from  this  model  are 
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where  m  and  e  are  the  electron  charge  and  mass,  and  g  is 
a  scaling  parameter  defined  in  [3].  In  deriving  these 
equations,  the  pressure  at  the  cathode  was  chosen  so  that 
the  previous  MITL  equations  with  Q  =  0  are  recovered 

[2-4], 

Typically,  Eq.  (6)  is  referred  to  as  the  pressure 
balance  equation  and  Eq.  (9)  is  referred  to  as  the  voltage 
equation.  When  Qc  =  0  (i.e.,  Ec  =  0),  Eqs.  (6)  and  (9)  can 
be  combined  to  yield  the  traditional  voltage  equation 
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which  is  used  to  calculate  the  voltage  from  measuring  Ia 
and  Ic.  Previously,  Eq.  (10)  was  derived  by  postulating  an 
expression  for  the  charge  density  in  the  flow  layer.  Here, 
it  is  derived  directly  from  the  fluid  model. 

The  model  also  yields 
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so  that  Z0  >  Zf  >  Z0/2  for  an  electron  layer  thickness  ds 
ranging  from  0  to  d.  Additionally,  the  charge  density  in 
the  electron  layer  p,  the  axial  fluid  velocity  in  the  electron 
layer  vz,  the  electron  layer  thickness  ds,  the  voltage  at  the 
edge  of  the  electron  layer  Vs,  and  the  magnetic  vector 
potential  at  the  edge  of  the  electron  layer  As  are  given  by 
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Finally,  Eqs.  (7)  and  (9)  can  be  combined  to  provide 
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which  for  cQc  =  0  again  yields  the  traditional  result  [3] 
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In  practice,  the  dynamic  MITL  model  will  receive 
advanced  values  of  V,  Ia,  and  A  along  the  line  for  each 
time  step  from  the  TLC,  where  A  is  advanced  using  Eq. 
(11  a)  from  [4].  Then,  Eqs.  (6)  -  (9)  are  used  to  solve  for 
cQa,  cQc,  Ic,  and  Zf  to  determine  the  flow  along  the  line  for 
the  subsequent  time  step.  The  continuity  equation 
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is  used  to  calculate  current  lost  to  the  anode,  where  Ja  is 
the  linear  current  density  and  Ibss  =  cxJa  is  the  loss  current 
for  a  TLC  element  of  length  x  (in  seconds). 


III.  SOLUTIONS 

The  TLC  model  for  MITL  flow  depends  on  a  robust 
numerical  technique  for  solving  Eqs.  (6)  -  (9).  The 
equations  can  be  normalized  to  the  local  vacuum 
impedance  Z0  of  the  particular  TLC  element  under 
consideration.  Thus,  Z0cQa,  Z0cQc,  Z0Ic,  and  Z/Z0  are  the 
unknowns  and  V,  Z0Ia,  and  cA  are  the  known  parameters. 
Equations  (6)  -  (9)  can  be  combined  to  give  a  quintic 
equation  in  Z0cQa,  which  may  have  1,  3  or  5  real  roots. 
Valid  roots  are  restricted  to  physically-meaningful  regions 
of  parameter  space. 

To  map  out  this  parameter  space,  one  may  focus  on 
solutions  in  which  Z0cQc  =  0  (from  traditional  equilibrium 


Figure  2.  Plots  of  Z0IC  and  cA  as  functions  of  Z0Ia  for  V  = 
6  MV  with  Z0cQc  =  0. 


Figure  3.  Plot  of  Z0Ic  solutions  as  a  function  of  Z0Ia  for 
the  case  with  V  =  6  MV  and  cA  =  7.5  MV. 

MITL  theory).  First,  the  traditional  Z0Ia  versus  Z0IC 
solutions  for  the  Z0cQc  =  0  is  displayed  in  Fig.  2  for  V  =  6 
MV.  Also  shown  is  cA  versus  Z0Ia,  which  is  derived  from 
Eq.  (8)  with  Z0cQc  =  0.  Although  both  Z0Ic  and  cA  are 
really  functions  of  Z0Ia,  the  axes  have  been  oriented  with 
the  independent  variable  on  the  vertical  axis  in  order  to 
display  the  results  in  the  traditional  manner.  There  is  a 
one  to  one  correspondence  between  these  two  solutions 
with  the  example  indicated  by  the  arrow  in  Fig.  2 
corresponding  to  the  solution  in  Fig.  3  for  cA  =7.5  MV. 
Over  most  of  the  range  of  cA,  only  one  Z0cQc  =  0  solution 
exists.  However,  over  a  narrow  range  of  cA,  there  is  a 
region  where  two  solutions  with  different  Z0Ia  values 
exists,  corresponding  to  the  region  left  of  the  minimum  on 
the  Z0Ia  versus  Z0Ic  curve.  This  is  traditionally  referred  to 
as  the  saturated-flow  region. 

As  an  example  of  general  solutions  (i.e.,  with  Z0cQc 
not  necessarily  equal  to  zero),  consider  solutions  for  V  =  6 


Z0cQc(MV) 

Figure  4.  Plot  of  Z0cQc  solutions  as  a  function  of  Z0Ia  for 
the  case  with  V  =  6  MV  and  cA  =  7.5  MV. 
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Figure  5.  Plot  of  Z0cQa  solutions  as  a  function  of  Z0Ia  for 
the  case  with  V  =  6  MV  and  cA  =  7.5  MV. 

MV  and  cA  =  7.5  MV  as  a  function  of  Z0Ia.  This  value  of 
cA  occurs  in  Fig.  2  in  a  region  where  there  is  only  oneZ0Ia 
solution  with  Z0cQc  =  0.  Solutions  for  Z0Ic  are  shown  in 
Fig.  3.  Again,  the  axes  are  oriented  to  show  the  familiar 
Z0Ia  versus  Z0Ic  curve.  The  straight  dashed  lines  indicate 
bounding  constraints  on  the  solutions;  thus,  only  solutions 
with  values  that  fall  within  the  triangular  region  are 
meaningful.  These  bounds  are  determined  by  Z0  >  Zf  > 
Z/2,  Ia  >  Ic,  and  the  constraints  derived  from  these 
conditions. 

Also  note  the  traditional  Z0cQc  =  0  solutions,  indicated 
by  the  dotted  curve  in  Fig.  3.  Imagining  that  there  is  a 
Z0cQc  axis  projecting  in  and  out  of  the  paper  and  that 
solutions  with  different  values  of  Z0cQc  are  projected  onto 
the  plane  of  the  paper,  the  only  solution  with  Z0cQc  =  0 
occurs  at  the  location  indicated  by  the  “x”.  Note  that  this 
solution  falls  on  the  dotted  Z0cQc  =  0  curve  and 
corresponds  to  the  solution  indicted  by  the  arrow  in  Fig.  2 
for  cA  =  7.5  MV.  Other  points  that  seem  to  cross  the 
Z0cQc  =  0  curve  actually  lie  off  the  plane  of  the  paper  with 
Z0cQc  ^  0. 
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Figure  6.  Plot  of  Z/Z0  solutions  as  a  function  of  Z0Ia  for 
the  case  with  V  =  6  MV  and  cA  =  7.5  MV. 


Solutions  for  Z0cQc,  Z0cQa,  and  Z/Z0  as  functions  of 
Z0Ia  are  shown  in  Figs.  4,  5,  and  6,  respectively.  For  ease 
of  comparison,  the  axes  are  also  oriented  in  the  same 
fashion  as  Fig.  2.  In  Fig.  4  the  only  Z0cQc  =  0  solution  for 
these  values  of  V  and  cA  occurs  on  the  lower  branch. 
Note  that  when  Z0cQc  >  0  (i.e.,  Ec  >  0),  electron  emission 
occurs  at  the  cathode,  and  when  Z0cQc  <  0  (i.e.,  Ec  <  0), 
electron  re -trapping  occurs.  Also  note  that  Z0Ic  and  Z0cQc 
change  rapidly  near  the  Z0cQc  =  0  solution  for  small 
changes  in  Z0Ia,  while  Z0cQa  and  Z/Z0  do  not.  This 
suggests  that  strong  emission  or  re-trapping  will  tend  to 
keep  solutions  in  the  vicinity  of  the  Z0cQc  =  0  solutions  at 
least  when  the  flow  is  on  this  lower  branch. 


IV.  SUMMARY 

A  fluid  model  has  been  developed  to  model  generalized 
MITL  flow.  By  explicitly  including  pressure  in  the 
model,  the  voltage  equation  can  be  derived  directly  rather 
than  resulting  from  postulating  an  expression  for  the 
density  in  the  flow  layer  as  was  previously  done.  By 
allowing  solutions  with  Ec  ^  0,  emission  and  re-trapping 
can  be  treated.  By  assuming  that  both  the  charge  density 
and  axial  fluid  velocity  are  uniform  in  the  flow  layer,  a 
well-defined  flow  impedance  is  derived  with  Zf  =  Zm. 
This  assumption  is  allowable  due  to  the  inclusion  of  a 
pressure  term  in  the  fluid  model.  And,  finally,  by 
choosing  the  form  of  the  pressure  term  appropriately,  the 
model  reproduces  the  new  rescaled  MITL  flow  model 
when  Ec  (or  equivalently  Qc)  is  taken  to  zero.  This  newly 
extended  MITL  flow  model  may  be  combined  with  time- 
dependent  field  equations  to  build  a  dynamic  model  for 
MITL  flow  in  a  TLC.  A  robust  numerical  technique  is 
being  constructed  to  solve  these  new  generalized  MITL 
flow  equations,  however,  techniques  to  treat  impedance 
transitions,  adders,  non-emitting  regions,  load  coupling, 
etc.  are  also  needed  to  complete  the  TLC  model. 
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